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Abstract. 

We investigate the effect of deviations from general relativity on approach to the 
initial singularity by finding exact cosmological solutions to a wide class of fourth-order 
gravity theories. We present new anisotropic vacuum solutions of modified Kasner type 
and demonstrate the extent to which they are valid in the presence of non-comoving 
perfect-fluid matter fields. The infinite series of Mixmaster oscillations seen in general 
relativity will not occur in these solutions, except in unphysical cases. 
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1. Introduction 

There have been numerous studies in the literature on generalisations of the usual 
Einstein-Hilbert action of general relativity to more complicated functions of the 
curvature (see refs. in PUE])- The motivation for such studies comes from a variety of 
different sources ranging from attempts to include quantum effects in the gravitational 
action [3] to the investigation of phenomena that are presently inadequately explained 
in the standard model, such as dark energy [HE], and violations of the cosmic no hair 
theorem jHllZj. Of particular interest is the behaviour of these modified theories of gravity 
in the high-curvature limit, where quantum corrections to general relativity are expected 
to become important and their influence on the occurrence of singularities and possible 
bounces of the universe at high curvatures can be studied jHj. It is the investigation 
of these modified theories of gravity on approach to the initial cosmological singularity 
that concerns us here. Previous studies have focussed on the behaviour of isotropic 
cosmologies but, as we know from the situation in general relativity, their behaviour can 
be misleading. Anisotropies diverge faster than isotropic densities at high curvatures and 
will generally dominate the behaviour of the cosmology at early (and in some cases even 
late) times. The most important anisotropic solutions of general relativistic cosmologies 
are the vacuum Kasner solutions [§| and their fluid-filled counterparts that form Type I of 
the Bianchi classification of three-dimensional homogeneous spaces. These universes are 
geometrically special, in vacuum (or perfect-fluid) cases they are defined by just one (or 
two) free constant (s) compared to the four (eight) that specify the most general spatially 
homogeneous vacuum (or perfect fluid) solutions. However, they have proved to provide 
an excellent dynamical description of the evolution of the most general models over finite 
time intervals. The chaotic vacuum Mixmaster universe of Bianchi Type IX undergoes 
a infinite sequence of chaotic space-time oscillations on approach to its initial or final 
singularities which is well approximated by a sequence of different Kasner epochs which 
form a Poincare return mapping for the chaotic dynamical system [TUllll[IT2llT51ITll ll5|. 

The way in which the Kasner metric has played a central role in the elucidation 
of the existence and structure of anisotropic cosmological models and their singularities 
in general relativity makes it the obvious starting point for an extension of that 
understanding to cosmological solutions of higher-order gravity theories. In an earlier 
Letter JE] we have reported the discovery of a new class of exact Kasner-like solutions 
for gravity theories derived from Lagrangians that are an arbitrary power of the scalar 
curvature. In this paper we will generalise that study to a much wider class of higher- 
order Lagrangian theories. We will determine the conditions for the existence of Kasner 
solutions and find their exact forms. In some cases these solutions are required to 
be isotropic and correspond to exact Friedmann-Robertson- Walker (FRW) vacuum 
solutions with zero spatial curvature. By studying gravity theories whose Lagrangians 
are derived from arbitrary powers of the curvature invariants we are able to find simple 
exact solutions. Past studies have usually focussed on the addition of higher-order 
curvature terms to the Einstein-Hilbert Lagrangian. This results in enormous algebraic 
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complexity and exact solutions cannot be found. The results presented here provide a 
tractable route into understanding the behaviours of anisotropic cosmological models 
in situations where the higher-order curvature corrections are expected to dominate. 
The solutions we present are vacuum solutions but we provide a simple analysis which 
determines when the introduction of perfect fluids with non-coming velocities has a 
negligible effect on the cosmological evolution at early times. 



2. Field Equations 

We consider field equations derived from an arbitrary analytic function of the three 
possible linear and quadratic contractions of the Riemann curvature tensor; R, R a bR ab 
and R a bcdR abcd The weight-zero Lagrangian density for such theories is 

C G = X ~ 1 V9f(X,Y,Z) (1) 

where x is a constant and f(X, Y, Z) is an arbitrary function of the variables X, Y and 
Z defined as X = R, Y = R a bR ab and Z = R abcd R abcd . We ignore the boundary terms, 
for a discussion of which see ref. [IB]. The field equations derived from the variation of 
the corresponding action are [TBI CHj 

P ab = - l -fg ah + f x R ab + 2f Y R ac R b c + 2f z R acde R b cde (2) 

- D(f Y R ab ) - g ab (f Y R cd )., cd + 2U Y R c{a ) b \ 

- fx-A9 ab 9 cd - 9 ac 9 M ) - 4(fzR d{ab)c ); Cd 

_ X rp 

- -lab 

where a subscript X, Y or Z denotes differentiation with respect to that quantity. 

We will be looking for spatially homogeneous, vacuum solutions of Bianchi type I, 
described by the line-element 

3 

ds 2 = -alt 2 + J2 t2Pt d* 2 i (3) 
i=i 

where p±, P2 and p% are constants to be determined. For the special case Pi = P2 = Pz 
these solutions correspond to spatially flat FRW metrics, which will be found to exist 
for various higher-order theories both in vacuum (where none exist in general relativity) 
and in the presence of a perfect fluid. 

We can determine the number of independently arbitrary functions of three space 
variables that will characterise the general vacuum solution of the field equations on a 
Cauchy surface of constant time in these higher-order theories. The field equations are 
in general 4 th order in time; so if we choose a synchronous reference system then we 
need 6 functions each for the symmetric 3x3 tensors g a p, g a p, g a p and g a p. This 
gives 24 functions, but they may be reduced to 20 by using the 4 constraint equations, 
and finally again to 16 by using the 4 coordinate covariances. If a perfect fluid were 
included as a matter source the final number would rise to 20 due to the inclusion of the 
density and 3 non-comoving velocity components in the initial data count. We note that 
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in general relativity the general vacuum solution is prescribed by 3 arbitrary functions 
of 3 space variables. In general relativity, the Kasner vacuum solution is prescribed by 
one free constant, as the three pi satisfy two algebraic constraints. 

3. Exact Solutions 

For a realistic theory we should expect the dominant term of the analytic function 
f(X,Y,Z) to be of the Einstein-Hilbert form in the Newtonian limit. However, there 
is no reason to expect such a term to dominate in the high curvature limit - in fact, 
this is the limit in which quantum corrections should become dominant. We therefore 
allow the dominant term in a power series expansion of f(X, Y, Z) to be of the form R n , 
(R a bR ab ) n or (R a bcdR abcd ) n on approach to the singularity, where n is a constant. These 
three different cases will be investigated separately below. 

3.1. f = R n 

In a previous work ^B] we found Bianchi type I solutions to the theory defined by the 
choice / = f{X) = R n . These solutions have a particularly simple form and they could 
be used to show that anisotropic universes in such theories do not exhibit an infinite 
sequence of Mixmaster oscillations on approach to the initial singularity, if n > 1. We 
also include these solutions here for completeness. 

Substituting / = f(X) = R n into the field equations (J2J) together with the line- 
element (J3J) gives the two independent algebraic constraints, 

(2n 2 - An + 3)P + (n - 2)Q - 3(n - l)(2n - 1) = 

and 

2(n 2 - 1)P + (n - 2)(P 2 - Q) - 3(n - l) 2 (2n - 1) = 0, 
where we have defined 

p = zZ Pi and Q = zZ p2 i- ( 4 ) 

i i 

The constraint equations have two classes of solution. The first is given by 
3(n- l)(2n- 1) 



P 



p2 

Q = — 

* 3 



(2-n) 



which is only solved by the isotropic solution p\ = pi = pz = P/3. This is the zero- 
curvature isotropic vacuum universe found by Bleyer and Schmidt [21 120] • It can be seen 
that these isotropic cosmologies are valid for all n ^ 2 and correspond to expanding 
universes for n<l/2orl<n<2 and to collapsing universes for l/2<n<lor2<r2. 
The second class of solutions to (JU) is given by 

P = 2n-\ (5) 
Q = (2n- 3)(1 - 2n). 
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n 




value of pi 



Figure 1. For the solution the intervals in which the Kasner indices pi are allowed 
to lie can be read off from this graph. For any value of n in the range | ^ n ^ | a 
horizontal line is drawn; the boundaries of the intervals in which the pi lie are then 
given by the four points at which the horizontal line crosses the two closed curves. For 
n = 1 these boundaries can be seen to be 0, | and 1, as expected for the Kasner 
solution of general relativity. 

This class of solutions represents a generalisation of the Kasner metric of general 
relativity and was the basis of the study It was shown in [TOj that this class 

of solutions only exists for n in the range 1/2 ^ n ^ 5/4 and that the values of the 
constants pi must then lie within the ranges 

2n - 1 - 2 v / (2n-l)(5-4n) 3p 1 ^ 2n - 1 - y/(2n- l)(5-4n), 
2n - 1 - y/(2n- l)(5-4ra) < 3p 2 ^ 2n - 1 + ^{2n - 1)(5 - An), 
2n - 1 + ^(2n - 1)(5 - An) < 3p 3 ^ 2n - 1 + 2^(2™ - l)(5-4n), 

where it has been assumed without loss of generality that p\ ^ p 2 ^ P3- These ranges 
are shown on Figure 1.111 and can be read off by drawing a horizontal line of constant 
n on the plot. The four points at which the horizontal line crosses the two closed 
curves gives the four boundary values for the allowed ranges of the constants pi. It can 
be seen that the points at which the curves cross the abscissa, which corresponds to 
n — 1, give the boundary values —1/3 ^ p\ ^ ^ P2 ^ 2/3 ^J53 ^ 1, in agreement 
with the Kasner solution of general relativity. For n > 1/2 these solutions correspond 
to expanding universes with a curvature singularity at t — 0. For n — 1/2 the only 
solution is Minkowski spacetime. 

For a universe filled with a perfect fluid having an equation of state p = (7 — l)p, 
7 constant, relating the fluid pressure p to its density p, the field equations (j2J) have the 
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isotropic solution 

Pl=P2=P3 = (6) 

for 7/fl. This reduces to the spatially-flat FRW solution of general relativity in 
the limit n — > 1. For n > 0, these isotropic cosmologies are expanding, and for 
n < they are contracting, with n = giving Minkowski spacetime. The stability 
and observational consequences of cosmologies of this type were investigated in |T§l l2"T]. 
where primordial nucleosynthesis and the microwave background were used to impose 
observational constraints on the admissible values of n. 



3.2. f= (R ab R ab ) n 

Substituting / = f(Y) = (R ab R ab ) n into the field equation (J2J) along with the metric 
ansatz (jHJ) gives the two independent equations, 

yn-l ^p2 + Q_ A pQ + p2 Q + g 2) _ (?) 

2(3P 2 + P 3 + 3Q - 9PQ + 2Q 2 )n + 2(4P 2 + 4Q - 8PQ)) = 

and 

Y n-i (j24P - 2P 2 - 2P 3 - 30Q + 10PQ)n + (8P 2 - 32P + 2AQ)n 2 (8) 

+ (P 2 - 4P + 2P 3 + 9Q- 10PQ - P 2 Q + 3Q 2 )) = 0, 

where P and Q are defined as before. Equations (J7J) and (jHJ) have five classes of solutions. 
The first of these classes is given by 

P = Q = 0. 

This class of solutions is only satisfied by p\ = P2 = P3 = 0, which is Minkowski 
spacetime. It can be seen from (JZJ) and (jEJ) that as this solution corresponds to Y = 
it only exists for n ^ (for n < the premultiplicative factor in ((7j) and (jHj) causes the 
left-hand side of those equations to diverge). 
The second class of solutions is given by 

P = Q=l. 

This is just the Kasner solution of general relativity, for which the values of the constants 
Pi are constrained to lie within the ranges — l/3^p\^0^p2^2/3^p^^l. Again, 
this solution corresponds to Y = and so is only valid for n ^ 0. 
The third class of solutions is given by 

_ 3(1 - 3n + An 2 ) ± ^(-1 + lOn - 5n 2 - 40n 3 + 48n 4 ) 
~ 2(1 -n) 




The only solution belonging to this class corresponds to an isotropic and spatially flat 
vacuum FRW cosmology. The values of the constants Pi are all equal to P/3 in this case 
and the solution is valid for all n 7^ 1. 
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Figure 2. For the solution the intervals in which the Kasner indices pi are 
constrained to lie can be read off this graph, as with figure A horizontal line 

of constant n is drawn; the boundaries of the intervals in which the pi lie are then 
given by the four points at which the horizontal line crosses the curves. For n = 1 
these boundaries can be seen to be — -|, 0, | and 1, as expected for the Kasner solution 
of general relativity. 



The fourth class of solutions is given by 

P = (1 - 2nf (9) 
Q = 1 - 8n + lQn 2 - 8n 3 . 

The solutions belonging to this class are anisotropic cosmologies which reduce to the 
standard Kasner form, P — Q — 1, in the limit n — > 1. For this class of solutions, the 
values of the constants pi are in general constrained to lie within the ranges 

(1 - In) 2 - 2A ^ 3pi < (1 - 2n) 2 - A 
(1 - 2n) 2 -A ^ 3p 2 ^ (1 - 2n) 2 + A 
(l-2n) 2 + A ^ 3p 3 ^ (1 - 2n) 2 + 2A 

where A = — 2n)(l — 6n + 4n 3 ) and the pi have again been ordered so that 
Pi ^ P2 ^ P3- These boundaries are shown in figure 13.21 which can be read in the 
same way as figure 13. 1| by taking a horizontal line of constant n and noting the four 
points at which that line intersects the curves. These intercepts determine the allowed 
intervals for the values of the p^. For real- valued Pi, the value of n must lie either in 
the range n\ ^ n ^ n 2 or in the range 1/2 ^ n ^ n s , where ri\, n 2 and n 3 are the 
roots of the cubic polynomial 1 — 6n + 4n 3 = and are chosen such that n\ < n 2 < n 3 . 
These generalisations of the Kasner metric always correspond to expanding cosmologies, 
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independent of the value of n. 

The fifth class of solutions to (JJJ) and (jHJ) is given by 

P = An - 1 

Q = -3 + 12n - 8n 2 ± 2(1 - 2n) ^2(1 - An + 2n 2 ). 

This class describes complex- valued pi for all values of n (except for n = 1/4, for which 
this class of solutions reduces to the first class) and is therefore of limited interest. 

The isotropic solution for a universe filled with a fluid with equation of state 
p = (7 — l)p is given by the choice 

An 

37 

This solution reduces to the spatially-flat FRW cosmology of general relativity in the 
limit n — > 1/2 and corresponds to an expanding universe for n > and to a contracting 
universe for n < 0. 



00 £ I T) T3abcd\n 

J.J. J — [Kabcd-tl ) 

Substituting / = f(Z) = (R abcd R abcd ) n and © into the field equations (0) gives the two 
independent equations 



n 



2 12 ^ 4 3 
+ (— + 6Q + 2PQ - 2P 2 Q + 3Q 2 - 105 + — ) 

+ 8(S-Q)rA =0 



and 

\v 4 222 2 2 4 / v; 

+ (6P + 4P 2 - 16Q - 2Pg + 85)n 

+ 8(Q - P)n 2 ^) = 
where P and Q are defined as before and we have also now defined 

i 

Equations (|10|) and ()11|) have four different classes of solution. 
The first class of solutions is given by 

P = Q = S = 0. 

The only solution that belongs to this class is Minkowski spacetime. As R a i )C dR abcd = 
for Minkowski spacetime it can only be a solution for n ^ 0, due to the premultiplier in 
equations (fTUj) and (fTTj). 
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The second class of solutions is given by 
P = Q = S=\. 

The only solution that belongs to this class corresponds to p\ = P2 = and p^ = 1. 
Making the coordinate transformations z = t sinh z and i = t cosh z [21] allows the 
line-element (J3J) to then be written in the form 

ds 2 = —dt 2 + dx 2 + dy 2 + dz 2 

which is clearly Minkowski spacetime again. This is only a solution for n > 0, as with 
the first class of solutions. 

The third class of solutions to (fTU|) and (fTTjl is given by 



p _ 3(1 - 2n + An 2 ± y/-l + 10ra - 16n 2 + 16n 4 
~ 2(1 - n) ' 

P 2 

^ = T' 

9 ' 

which only has the isotropic and spatially flat vacuum FRW cosmology as a solution, 
where pi = P2 = P3 = P/3. This is a solution for all n ^ 1. 
The fourth and last class of solutions is given by 

P = 4n-l 

Q = -{16n 2 - 8n - 1 ± 4^2^(1 - 2n) + S(l -n)]}. 
3 

This class of solutions corresponds to anisotropic metrics with Z = and is unusual in 
that it cannot be expressed in the form P = constant and Q = constant. This feature 
means that the standard picture of a plane intersecting an ellipsoid is no longer a valid 
one for this class. Solutions in this class do not appear to have any range of n for which 
the constants pt take real values, and so are of limited physical interest. 

The isotropic solution for a universe filled with a fluid with equation of state 
p = (7 — l)p is given by 

4n 

Pl=P2=P3 = ir- 

As before, this reduces to the spatially flat FRW solution of general relativity in the 
limit n — > 1/2 and corresponds to an expanding universe for n > and to a contracting 
universe for n < 0. 



4. Investigation of the effects of matter 

We have identified above two generalisations of the Kasner metric of general relativity 
which are solutions to the scale-invariant theories of gravity we are investigating. The 
first of these is a solution to the theory defined by C = R n and was the subject of 
investigation in [TH]; the second is a new solution to the theory C = (R a bR ab ) n . These 



Further Exact Cosmological Solutions to Higher-Order Gravity Theories. 



10 



are the first exact anisotropic solutions to be found for higher-order gravity theories. In 
this section we will investigate some of the properties of these cosmologies. 

The behaviour of these solutions is of particular interest when considering the 
Bianchi type VIII or type IX 'Mixmaster' cosmologies. The field equations for these 
cosmologies can be cast into the form of the equations of motion of a particle moving 
inside an exponentially steep triangular potential well with open channels in the corners 
[TT] . The three steep-sided walls are created by the 3-curvature anisotropies. In the 
region where the potential is negligible (far from the walls) the behaviour of the solutions 
approaches that of the Kasner metric. As the exponentially steep potential wall is 
approached the universe 'particle' is reflected and re-enters a Kasner-like regime with 
the Kasner indices Pi systematically changed to some new values by the rule of reflection 
from the potential wall. In general relativistic cosmologies of Bianchi types VIII and 
IX this process is repeated an infinite number of times as the singularity is approached 
[HH [TTJ H2 HS1 UH EES] so long as matter obeys an equation of state with p < p. After 
reflection from the potential wall the Kasner index that was previously negative is 
permuted to a positive value and the lowest-valued positive Kasner index is permuted 
to a negative value. This is repeated ad infinitum in the general relativistic cosmology 
as one of the Kasner indices must be kept negative while the other two are positive. In 
the generalisations of the Kasner metric presented above this is no longer the case; for 
some solutions it is possible for all of the Kasner indices to be positive. If the Kasner 
indices are being permuted in a chaotic fashion for long enough then, eventually, they 
will end up in such a configuration. Once this occurs the oscillations will end, as all 
spatial dimensions will be contracting, and the singularity will reached monotonically 
without further oscillation of the scale factors [22| 123]. 

From figure 13.11 it can be seen that for solution (jSJ) all Kasner indices can be made 
positive for values of 1 < n < 5/4, in the C = R n theories. Similarly, from figure l3~2l 
it can be seen that for the solution (jH)), with n\ < n < n2 or 1 < n < n%, that all the 
Kasner indices can be made positive, in the theory C = (R a bR ab ) n (where the are 
defined as before). For both of these solutions, with these ranges of n, there is therefore 
no infinite series of chaotic Mixmaster oscillations as the singularity is approached. 
However, the solutions to both of these theories are still expected to exhibit an infinite 
number of chaotic oscillations if 1/2 < n < 1, as at least one of the Kasner indices 
must be kept negative whilst another is kept positive. A separate detailed analysis is 
require to determine where there are differences between this behaviour and the chaotic 
oscillations found in general relativity. 

It now remains to show that the vacuum solutions found above are the asymptotic 
attractor solutions in the presence of non-comoving matter motions as t — > 0. This 
analysis will follow closely that of ^H] and [21]. We aim to show that the fluid stresses 
diverge more slowly than the vacuum terms as t — > and so produce negligible metric 
perturbations to an anisotropic Kasner universe. If this is the case, then the vacuum 
solutions above can indeed be considered as the asymptotic attractor solutions even in 
the presence of matter. Matter will just be carried along by the expansion and behave 
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like a test fluid. 

We now consider a perfect fluid with equation of state p = (7 — l)p where 1 ^ 7 < 2 
and non-zero 4- velocity components, u 1 normalised so that u l Ui = 1. The conservation 
equations T ab .f, = on the metric background J3J) can then be written in the form 

jpu k (u iik - -u l g M ,i) = - (7 - l)uiU k p )k . 

Neglecting spatial derivatives with respect to time derivatives these equations integrate 
to 

1 

t p uop~> = const. 

(7-1) 

u a p -1 = const. 

The neglecting of spatial derivatives means that these equations are valid on scales larger 
than the particle horizon as in the velocity-dominated approximation (although we have 
not so far restricted the fluid motions to be non-relativistic) |2H ESI EZ| • 

From the second of the integrals above, it can be seen that all the covariant 
components of the spatial 3-velocity, approximately equal. This is not true of 

the contravariant components as the Kasner indices, Pi, in the metric elements used to 
raise indices are not equal in these solutions. The contravariant component that diverges 
the fastest, and dominates the others in the t — > limit, is therefore u 3 = ust 2m , as 
p 3 is the largest of the pj. If the 4- velocity is normalised, so that u a u a = 1, and the 
contravariant 3-velocity component u 3 diverges the fastest as t — > 0, then we must have 
UqU° ~ M3M 3 = (u 3 ) 2 t~ 2m in that limit. The integrated conservation equations above 
can now be solved approximately in this limit to give 

p ^ £-7(Pi+P2)/(2-7) 
Ua ^ £(pi+pa)(7-l)/(2-7) 

as t — > 0. 

It is now possible to calculate the leading-order contributions to the energy- 
momentum tensor, T 6 a = (p + p)u a Ub — pS^, as being 

T ° ~ P u°u ~ r p - pi 

Tl ~ p ~ r 7(P-P3)/(2- 7 ) 

r| ~ P u 2 u 2 ~ r 2p2 - p+P3 
t| ~p M 3 M 3 ~r p - p3 . 

The component which diverges the fastest here as t — > is T 3 3 ~ t~ p ~ P3 , for general 
7. For the case 7 = | where the 4- velocity of the fluid is comoving, «j = <5 4 °, all the 
components of T fe a diverge as p ~ f- 2 P+ 2 P3_ 

For the generalisation of the Kasner metric in the £ = R n theory, (jSJ), the vacuum 
terms diverge as t~ 2n . We therefore require that 2n > P + p 3 , or p 3 < 1, in order 
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for the vacuum terms to dominate over the matter terms for a non-comoving perfect 
fluid, with general 7, in the limit t — > 0. For the 7 = |, comoving perfect fluid 
(ui = 5®), the condition for vacuum domination is 2n > 2P — 2p%, or p% > n — 1. 
Both of these conditions are ensured by the boundary values on £>3 given above. In the 
velocity-dominated limit, the vacuum solutions given here are therefore the appropriate 
asymptotic solutions on approach to the singularity, as was found in [IE] . 

For the generalised Kasner metric in the L = {R a i 1 R ab ) n theory, Q, the vacuum 
terms diverge as t~ 4n . For vacuum domination as t — > in this solution we therefore 
require 4n > P + p 3 , or p 3 < — 1 + 8n — 4n 2 , for the general non-comoving 7 fluid; and 
An > 2P — 2p%, or p 3 > 1 — Qn + 4n 2 , for the 7 = 4/3 comoving fluid. From figure PP1 
it can be seen that solutions of this class can be in one of two regions, n\ ^ n ^ n 2 
or 1/2 ^ n ^ n%. The validity of the vacuum solution as t — > is different in these 
two different regions. For the region 1/2 ^ n ^ n% the boundary conditions on the 
index p 3 mean that the conditions above are automatically satisfied. For the region 
n± ^ n ^ n 2 these inequalities are not always satisfied. For general 7, the condition for 
vacuum domination is only met if n lies in the narrow range 1/6 < n < n 2 , for any other 
value of n, in this region, the index p% can be such that the fluid diverges faster than 
the vacuum terms. Similarly, for the 7 = 4/3 comoving fluid the condition for vacuum 
domination is only satisfied if n lies in the narrow range n^ < n < n 2 , where n^ is the 
real root of the cubic 

4n 3 - 12n 2 + lOn - 1 = 0. 

To understand the evolution of an anisotropic solution of the form (J3J) in the theory 
C = (R a bR ab ) n , where n\ ^ n ^ n 2 , it is therefore necessary to take into account the 
relativistic motions of any fluid that is present (except in the narrow ranges of n identified 
above). This range of n can, however, be regarded as not belonging to physically 
interesting theories on other grounds. In order to agree with weak-field experiments it is 
necessary for a gravitational theory to contain a term which approximates the Einstein- 
Hilbert action in the weak-field limit (see ref. ^H])- Any theory containing such a limit 
must therefore have a term in its Lagrangian that diverges as t~ 2 . In considering the 
behaviour of alternative theories close to the singularity it is necessary for the extra 
terms in the Lagrangian to diverge faster than this if they are to be influential as t — > 0. 
For the theory C = (R a bR ab ) n this requires n > 1/2. For n < 1/2 such a term will 
diverge slower than the Einstein- Hilbert term which will then dominate and display the 
standard Kasner behaviour of general relativity. For the anisotropic solution Q we 
are therefore only interested in the region lying in the range 1/2 ^ n ^ 713, which has 
vacuum terms that dominate the fluid stresses in the velocity-dominated approximation. 

It remains to investigate the effects of fluids with stiff equations of state, 7 > 4/3. 
For sufficiently stiff fluids the velocity-dominated approximation is not valid and 
u a u a — > as t — > due to the very high inertia of the fluid producing a slow down 
under contraction [2Bj- In such a limit uq — > 1 and the conservation equations can be 
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solved to give 

p ~ r 7P 

In this approximation we can write « a M a ~ (u^) 2 t^ 2p3 ~ £ 2 (7-i)-P-2p 3 _ can now ^ e 
seen that this behaviour occurs when 

p 3 < (7 - (12) 

In this limit it is only required that p diverges more slowly than the vacuum terms, as 
pu a u a « p when t — > 0. 

For the solution (jSJ) to the theory £ = R n it can be seen from the condition (|12|) 
that it is necessary for 7 ^ 4/3 in order for the velocity-dominated approximation to 
break down (this is derived using the upper limit on n and the lower limit onp 3 ). The 
condition that the fluid stresses diverge more slowly than the vacuum terms is now 

2n . . 

1 < (2^T) < 13 > 

where it has been assumed n > 1/2. A sufficient condition to satisfy this for all allowed 
values of n is 7 < 5/3. For fluids with a stiffer equation of state, 5/3 ^ 7 < 2, the 
vacuum terms dominate over the fluid motions provided that n < 1. For n < 1, however, 
the Einstein-Hilbert term will be the leading one in the gravitational action, as discussed 
above. 

Similarly, for the solution Q to the theory C = (R a bR ab ) n it can be seen from 
the condition (jl2)l that it is again necessary for 7 ^ 4/3 in order for the velocity- 
dominated approximation to break down. In the region 1/2 ^ n ^ n 3 the condition for 
the vacuum solution to be unperturbed by the stiff fluid is always satisfied for any fluid 
with equation of state 7 < 2. However, in the region n\ < n < the condition for the 
vacuum to diverge faster than the energy density of the stiff fluid is never satisfied for a 
fluid which satisfies the necessary condition for the velocity-dominated approximation 
to break down (except in the narrow range (7 — v / 33)/8 < n < n 2 where there are some 
values of 7 for which the vacuum terms can diverge fastest). 



5. Conclusions 



We have investigated some anisotropic cosmological solutions to higher-order Lagrangian 
theories of gravity. Whist the standard general relativistic theory, defined by the 
Einstein-Hilbert action, appears to be consistent with all weak-field tests, there is 
less reason to think that it should be valid in high curvature regimes, such as in the 
vicinity of a possible initial cosmological curvature singularity. In fact, it is in the high- 
curvature limit that quantum effects should become important and we should expect 
to see deviations from the standard theory. Without knowing the exact form of such 
deviations, we have approached this problem by considering a general class of theories 
that can be derived from an arbitrary analytic function the three curvature invariants 
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R, R a bR ab an d R a bcdR abcd ■ Expanding this function as a power series in these variables 
we then expect the dominant term in the Lagrangian to be of the form R n , (R a f 1 R ab ) n 
or (R a b c dR abcd ) n as the singularity is approached. We have found all of the solutions 
to these theories that can be expressed in terms of the Bianchi type I line-element 
(JHJ) in vacuum. These solutions provide simple testing grounds for the exploration of 
quantum cosmological effects in higher-order gravity theory. We have also found the 
homogeneous, isotropic and spatially flat solutions to these theories in the presence of 
a perfect fluid. 

Exact vacuum anisotropic solutions of the form Q were found for the theories 
£ = R n and £ = (R a bR ab ) n , whilst it was shown that no such solutions exist for theories 
defined by £ = (R abcd R abcd ) n . The properties of these solutions, with respect to their 
relation to the more general Bianchi type VIII and IX cosmological behaviour, has been 
investigated. We have argued that for all of the physically relevant new solutions of these 
theories, the universe will not experience an infinite number of Mixmaster oscillations 
as the singularity is approached. This is an extension of what was found in the earlier 
work [T^]. The extent to which these vacuum solutions can be considered as realistic 
in the presence of a perfect fluid has also been investigated. It has been shown in the 
velocity-dominated approximation that all the anisotropic vacuum solutions found for 
plausible theories are valid in the limit t — > 0. The case of stiff fluids that do not 
satisfy the velocity-dominated approximation as t — > has also been investigated. For 
the solutions to the theories £ = R n it was found that for a fluid with equation of 
state 7 < 5/3 the vacuum solutions are good approximations in the vicinity of the 
singularity. For the theories £ = (R a bR ab ) n it was found that the vacuum solutions are 
good approximations for all 7 < 2. 
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